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Abstract. Spatial entropy redistribution plays a key role in the adiabatic 
cooling of ultra-cold lattice gases. We show that high-spin fermions with a 
spatially variable quadratic Zeeman coupling may allow for the creation of an 
inner spin-1/2 core surrounded by high-spin wings. The latter are always more 
entropic than the core at high temperatures, and remarkably at all temperatures 
in the presence of frustration. Combining thermodynamic Bethe Ansatz with local 
density approximation, we study the spatial entropy distribution for the particular 
case of one-dimensional spin-3/2 lattice fermions in the Mott- phase. Interestingly, 
this spatially-dependent entropy opens a possible path for an adiabatic cooling 
technique which contrary to previous proposals, would specifically target the spin 
degree of freedom. We discuss a possible realization of this adiabatic cooling, 
which may allow for a high-efficient entropy decrease in the spin-1/2 core and help 
accessing antiferromagnetic order in experiments on ultracold spinor fermions. 
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1. Introduction 

Ultra-cold atoms in optical lattices offer an extraordinary controllable scenario for 
the study of strongly-correlated systems [TJ [2] , exemplified by the observation of the 
superfhiid to Mott-insulator (MI) transition in ultra-cold bosons [3J. Remarkable 
progress has been achieved in lattice fermions as well, allowing for the precise 
analysis of the Fermi-Hubbard model, a key model in condensed-matter physics of 
particular relevance in the study of high-temperature superconductivity [3] • Exciting 
recent experiments have reported the realization of the metal to MI transition in 
two-component fermions [5J [B]. Due to super-exchange, the MI phase of spin- 1/2 
fermions is expected to exhibit a magnetic Neel (antiferromagnetic) ordering. On- 
going experiments are already very close to reach the mean- field entropy per particle, s, 
for Neel ordering in a 3D cubic lattice (s/fcs = In 2) [7]. However quantum corrections 
reduce the critical s down to s^/fes ~ 0.35 [HE]. 

Reaching such an extraordinary low entropy constitutes nowadays a major 
challenge, which demands novel types of cooling especially designed for many-body 
systems in optical lattices [10]. A number of cooling proposals have been recently 
suggested [TTJ El HH H3I HH H3 HSl H3, most of them based on the redistribution 
of entropy within the trap, where certain regions act as entropy absorbers from the 
region of interest, i.e. a Mott insulator at the trap center. 

Interestingly, spin degrees of freedom may be employed for designing cooling 
techniques resembling adiabatic demagnetization cooling in solid-state physics |18j . 
In this method, a decrease in the strength of an externally applied magnetic field 
allows the magnetic domains of a given material to become disoriented. If the 
material is isolated, temperature drops as the disordered domains absorb thermal 
energy in order to perform their reorientation. In cold atoms, this technique was 
pioneered in Chromium experiments, where the spin-flip mechanism was provided by 
dipolc-dipole interactions |19j . Recently, a novel demagnetization cooling mechanism 
has been proposed for two-component fermions based on time-varying magnetic field 
gradients [20]. In that method, scalar domains are cooled by transferring particle-hole 
entropy into magnetic entropy in overlapping regions between the two components. 
Note, however, that gradient cooling does not address cooling of the spin degrees of 
freedom, contrary to the method discussed below. 

In this paper, we study the spatial entropy distribution of multi-component spin- 
s' fermions with an inhomogeneous quadratic Zeeman effect (QZE), and how this 
spatially-dependent entropy profile may be employed for designing an adiabatic cooling 
method which specifically targets the reduction of spin entropy. We show in particular, 
that an inhomogeneous QZE may lead to an effective pseudo-spin-1/2 core surrounded 
by spin-S" fermions at the wings. We show that, remarkably, the spin-S 1 wings act as 
entropy absorbers all the way to vanishing temperatures in the presence of frustration. 
We illustrate the idea with the specific example of one-dimesnional (ID) spin-3/2 
lattice fermions in the Mott phase. Interestingly, the adiabatic growth of the spatially 
dependent QZE combined with spin redistribution via spin-changing collisions may 
open promising perspectives towards adiabatic cooling of the spin-1/2 core, which 
may in this way enter the antiferromagnetic spin coherent regime. 
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2. Entropy and frustration 

In the following we consider multi-component fermions loaded in an optical lattice. 
Increasing the number of spin components from 2 to TV (effective spin S = (N — l)/2) 
increases the capacity to allocate larger entropy per spin at high T by a factor 
of lnJV/ln2. This guarantees that at high T, well over the Neel ordering, spatial 
regions with a larger effective S act as entropy absorbers in the cooling process shown 
below. However, this simple argument does not apply at low T if the system acquires 
conventional Neel ordering [35]. The entropy of the Heisenberg antiferromagnet (HAF) 
in d dimensions scales for T -c T N (Neel temperature [2"5]) ass~ S~ d , and hence the 
entropy of a spin-S* system decreases as compared to that of a spin-1/2 HAF. This 
is clear since larger spins attach more to the Neel direction, and hence the density of 
states is smaller at low T leading to a lower entropy. 

The situation is reversed in frustrated systems, which present a large degeneracy 
of classical ground states with many branches of soft excitations at low T . One arrives 
at the simplest frustrated large- S model starting from the SU (N) symmetric Hubbard 
model 

H = -t ( c l, t c m,j+h.c.) + U/2j2nl (1) 

where m = (— S, . . . , S), c m ^ annihilates fermions with spin m in the site i, rij = 
^ m n m j = J2 m c m i c m.i, and t and U are the hopping and interaction coupling 
constants, respectively. In the strong-coupling limit, U 3> \t\, and retaining one 
fermion per site, one can derive the effective permutation model in second order of 
perturbation theory, 

H = J/2 p i,3> (2) 
<i,j> 

with Pij being the permutation operator and J — At 2 /U . This model is the SU(N) 
generalization of the HAF. It is exactly solvable in ID |27j and has N — 1 gapless spin 
modes, each dispersing at low momenta with velocity 2vin/N, where V\/2 — ft J/2 is 
the spin wave velocity of the spin 1/2 HAF. At low T, the entropy SU(N) spin model 
is larger than that of spin-1/2 HAF, due to increasing fluctuations by 'orbital' degrees 
of freedom, following s s (T) = N(N - 1)ttT/6v 1/2 [H]. For equivalent 2D and 3D 
lattice models an increase of entropy with unbinding number of degrees of freedom is 
also expected at low T. There, the classical ground state shows extensive degeneracy, 
and it is believed [2S] that the Neel order does not get stabilized for SU(N > 2) due 
to high frustration [2"5] . 

Hence, due to frustration, with increasing S, Neel order gets suppressed leading 
to a larger entropy storage capacity at low T. Thus, if the Mott edges get frustrated 
while preserving 5 = 1/2 character in the central region, one can use the frustrated 
edges as entropy absorbers all the way from high to extremely low T. 

3. ID Spin-3/2 fermions in the Mott-phase 

In the following, we illustrate the possibilities provided by high-spin lattice fermions 
with the specific case of a balanced mixture of spin-3/2 fermions in a ID optical lattice, 
in which the number of particles N m with spin m satisfies N m — N- m . Interparticlc 
interactions are characterized by the s-wave scattering lengths for channels with 
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Figure 1. Ratio between the entropy per particle so a * <? = (continuous line) 
or at q = q(T) (dashed line) (q(T) is the QZE at which entropy is largest at a 
given T), and the entropy s c at high QZE (q = qo), versus the temperature T in 
the lattice. Observe the crucial role of frustration at low T. If we were dealing 
with an unfrustrated 5 = 3/2 HAF, the curve would have arrived at low T at 
~ 0.4 [26] instead of 6. 

total spin and 2, ao,2- For s-wave interacting fermions S — 3/2 is the lowest 
spin allowing for spin-changing collisions (which preserve the total magnetization 
but transfers atoms from ±1/2 into ±3/2 and vice versa). Due to the conserved 
magnetization the linear Zeeman effect does not play any role. However, the QZE, 
characterized by the externally controllable constant g, induces a finite chirality 
t = t[(N 3/2 + N_ 3/2 ) - (JV 1/2 + N_ 1/2 )}. 

For large-enough interactions and at quarter filling (one fermion per site) the ID 
system enters into the Mott insulator regime, for which the ground state properties 
under QZE were studied in Ref. [31]. For large QZE the ground state is a pseudo spin- 
1/2 isotropic HAF. Reducing q < q cr (q cr = J In 2/2 at the ao = a 2 ) the system enters 
either a spin liquid phase (for a 2 < ao) or a dimerized phase (for a 2 > ao). For ao — a 2 
(the typical situation unless ao,2 are externally modified), the gap of the dimerized 
phase is exceedingly small, and hence the system behaves in practice as a spin liquid 
down to extremely low temperatures. For ao = a 2 , in the presence of a spatially 
variable QZE, the model Hamiltonian becomes H = H Q + [ J -m,i n m.i, with H given 
by Eq. © with N — 4. For homogeneous fi m ^, this model is exactly solvable and 
its thermodynamic properties may be calculated by means of thermodynamic Bcthc 
Ansatz. We follow the method of Ref. [32], based on the self-consistent solution of 14 
coupled integral equations, to obtain the corresponding free energy /. The chemical 
potentials for each component are li m ,i — M + m2 9ii where /i is the global chemical 
potential, and qi denotes the QZE constant at site i. The entropy is then given by 
s = —df/dT and the chirality by r = —df/dq. 

As mentioned above, for q > q cr the system becomes pseudo-spin-1/2. Hence, the 
ratio between the entropy per spin for q = (so) and that for q > q cr should follow 
the same dependence as s 3 / 2 (T)/ Si/ 2 (T). We illustrate this point in Fig.Q] where we 
compare s at large q = qo = 5J (s c ) to so- Note that at large T so/s c — 2, whereas at 
low T frustration leads to so/s c = 6. 
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Figure 2. Entropy per particle (in ks units) and chirality profiles for L = 
120 particles for a Gaussian QZE profile q(x) = qo exp (— 100x 2 /L 2 ) , with 
qo = 5J. Dashed lines indicate the initial entropy prior to the switching of the 
inhomogeneous QZE profile. Before the lattice loading, Ti/Tp =0.1 (left) and 
0.016 (right). Note that for q > q cr , t ~ 1 and the system retains a spin-1/2 
character. Note also how the entropy gain in the Mott core is larger for small T. 
The entropy bump for small T at q = q cr reflects the Van-Hove singularity at the 
bottom of the depleted Hubbard band. 

4. Spatially-dependent QZE 

We consider at this point the entropy and chirality profiles for the case of a 
non- homogeneous QZE, which may be achieved by means of microwave or optical 
techniques [22j [23l [24]. We perform local QZE approximation (similar to the local 
density approximation standard in trapped gases), i.e. we solve for the free energy 
at different positions x by varying q(x). The local QZE approximation demands a 
sufficiently slow variation of the QZE at the scale of the inter-site spacing. In this way 
we can evaluate the entropy profile inside the Mott insulator region. Note finally that 
although the calculation is done for the case ao = a,^ , the conclusions may be extended 
to the actual case in which ao and ai are slightly different, and spin redistribution via 
spin-changing collisions occurs. 

Figure [5] shows chirality and entropy profiles. Note that the entropy per site (i.e. 
per particle) is significantly larger at the Mott wings than at the center. Hence, if 
the total entropy is conserved in a process in which an initially spatially-independent 
entropy (deshed line on Fig. [2]) is brought to the profile of Fig. [2j the outer regions 
will remove entropy from the central core. This process is even more effective at low 
T . One may indeed estimate the entropy reduction (for T — > 0) at the Mott center 
for model ([2]) in ID, when considering an initial uniform spin-1/2 system and a final 
step-like distribution with a spin-1/2 core and spin S = (N — l)/2 wings: 

7 = Sl /s c = [(1 - L /L)N(N - l)/2 + L /L], (3) 

where Lq/L is the ratio of the number of sites in the spin-1/2 core to the total number 
of sites, and s, (s c ) is the initial (final) entropy per particle at the center. The gain is 
hence much bigger than the one at T > Tjy, which is obtained after one substitutes 
N(N - l)/2 by In N/ In 2 in Eq. ©. 
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Figure 3. Sliding down the isentropic curves (solid curves) with reducing 
uniformly QZE induces the reduction of the temperature in the lattice at the 
second stage of the proposed adiabatic cooling scheme. 

5. Adiabatic cooling 

The spatially-dependent entropy distribution opens interesting perspectives for 
adiabatic spin cooling. A possible scheme would consist on three steps. On a first 
stage a two-component balanced mixture is created at the lowest possible temperature 
Ti, being stabilized against spin-changing collisions |21j by means of a sufficiently 
large homogeneous QZE. On a second stage a lattice is adiabatically grown and 
the homogeneous magnetic field is adiabatically decreased allowing for spin-changing 
collisions throughout the sample. The drop of the temperature in the lattice with 
the adiabatic decrease of the homogeneous magnetic field can be estimated from 
the isentopic curves on Figure El However, to enter the spin coherent regime local 
entropy should be reduced. This is achieved in the final step which consists on 
slowly changing the QZE into a non-uniform profile by means of microwave or 
optical techniques [231 H31 121], leading to the coexistence of a spin- 1/2 Heisenberg 
antifcrromagnet (HAF) at the trap center and a spin- 5 spin liquid at the wings (Fig. [2]). 
Figure [4] shows, as a function of the initial Ti/Tp (Tp is the Fermi temperature 
of the original spin- 1/2 prior to the lattice loading, see below), the central entropy 
per particle with a high homogeneous QZE (dashed line) and with a Gaussian QZE 
(solid line). For the case considered s/fcs = 0.35 is achieved for the homogeneous 
case at Ti/Tp ~ 0.035 whereas for the inhomogeneous QZE profile it is reached at 
Ti/Tp ~ 0.09, showing that an inhomogeneous QZE may allow for a large entropy 
reduction at the center as a result of the entropy excess at the storage wings. Note 
that for Ti/Tp ~ 0.1 the central entropy may be reduced by more than factor of 2 
with respect to the entropy expected for the homogeneous spin- 1/2 case. 

Let us briefly comment on higher dimensional (d = 2,3) systems. A simple 
estimate may be obtained from a spin- wave analysis when ai > oq, where at least in 
the limit 3> ao Neel order sets in on bipartite lattices |33j , however due to frustration 




Figure 4. Entropy at the center s c (in ks units) as a function of the temperature 
before lattice loading, Ti/Tp, for the inhomogeneous q(x) profile of Fig. [T] (solid) 
and for a homogeneous one q = qo (dashed), for which the whole system retains 
a spin-1/2 character. In order to associate a central s c with a given Ti/Tp, 
we calculated for different T inside the lattice the entropy profiles and the total 
entropy Stot/L = ir 2 k B Ti/Tp. 



the modulus of the Neel order parameter decreases with increasing spin. There is just 
one spin- wave mode at q — > oo and three spin- wave modes at q = 0. Taking the values 
of spin- wave velocities from [33], we obtain that at T — s- the factor 7 is given by 
Eq. fl3j, albeit with N(N -l)/2 changed into 1 + 2[(a 2 + a )/(a 2 - a )] d , implying that 
the cooling should be even more efficient for higher dimensions. Note that the spin- 
wave analysis predicts 7 to increase indefinitely when approaching a 2 = oq, although 
spin- wave analysis becomes less reliable in the vicinity of that point [34] . 

6. Experimental feasibility 

A possible experimental sequence may be devised for e.g. 40 K. Prior to the lattice 
loading, a balanced mixture of e.g. F = 9/2, Mp = -5/2 and F = 9/2, Alp = -7/2 
states is prepared at the lowest possible Ti/Tp using standard techniques. At this 
stage a sufficiently large homogeneous magnetic field guarantees that the initial 
two-component mixture is stable against spin-changing collisions. The gas can be 
approximated with a good accuracy by a free Fermi gas and therefore the initial 
entropy per particle is given by s — kg-n Ti/Tp [3D]. State of the art experiments may 
reach at this stage Ti/Tp « 0.1, corresponding to an entropy per particle of s s=s kg- 
Once the gas is cooled down, a 3D lattice is grown and under proper conditions a Mott 
insulator with one particle per site develops at the trap center [5] [B] . The next step 
consists in slowly lowering the magnetic field to allow for quasi-resonant spin-changing 
collisions [35] throughout the sample leading to a redistributed population (between 
F = 9/2, Mp = -3/2, -5/2, -7/2, and -9/2 [36 ) and a significant drop of T in 
the lattice. Temperature drop after this stage for ID case, (assuming adiabaticity) 
is depicted on Fig. 3. Note that fermions with S > 1/2 suffer of large three-body 
loses [37] • However, in our scheme S > 1/2 is explored only at later stage, when the 
system is already in the hard-core regime with one fermion per site, thus three-body 
loses are largely suppressed. 

Finally, the F = 9/2, Mp = -3/2 and F = 9/2, Mp = -9/2 states may be slowly 
expelled from the trap center by the use of two Raman beam pairs. These beams 
couple F = 9/2, Mp = -3/2 to F = 7/2, Mp = -1/2 and F = 9/2, Mp = -9/2 
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to F = 7/2, Mp = —7/2. This AMp = — 1 coupling may be realized with just three 
lasers. A slight blue detuning ensures that the resulting dressed states experience a 
repulsive potential. The effect on the Mp = —7/2, —5/2 states is largely suppressed 
by the ratio between Raman detuning and Zeeman splitting, which is kept on the 
order of q cr . Furthermore, the Raman beams shift the resonance condition for spin- 
changing collisions in the center of the trap such that no more atoms in Mp = —9/2 
and Mp = —3/2 are produced. The outer regions have then a larger effective spin 
and hence act as entropy absorbers, as discussed above. Note that, interestingly, the 
center of the Raman beams acts as a crystallization point for a slowly-growing Neel 
state (if sjv is reached). In this way the problem of approaching the ground state 
despite the presence of many metastable low energy states may be circumvented. 

Let us mention some final remarks. For simplicity of our calculations we have 
considered cooling only within a Mott insulator. In the presence of particle-hole 
excitations we expect an even larger cooling efficiency [B]. Note that the S = 1/2 
region is deeper in the Mott region than the S > 1/2 region, and hence the S > 1/2 
region is expected to contain a larger particle/hole entropy. We stress also that the 
cooling efficiency is based on the assumption of adiabaticity. Spin-changing collisions 
are crucial in our cooling scheme. Although they are low energetic, they will be 
typically faster than super-exchange, and hence the adiabatic requirements of our 
model are comparable to those in other cooling schemes based on entropy relocations 
in optical lattices. 

7. Conclusions 

In summary, we studied a possible route for adiabatic spin cooling using high-spin 
lattice fermions. The process resembles demagnetization cooling, but the role of the 
magnetic field is played by a spatially dependent QZE, and spin flip is substituted by 
spin-changing collisions. The spatially dependent QZE leads to two distinct regions of 
different effective spin (S — 1/2 and S > 1/2). At high T the outer spin-S 1 region acts 
as an entropy absorber simply due to the larger spin, whereas the same remains true at 
even very low T due to frustration. As a result we showed that a significant reduction 
of the entropy (more pronounced at lower T and higher dimensions) of the spin-1/2 
Mott central region can take place. Magnetic refrigeration using spatially variable 
QZE can hence significantly facilitate experimental realization of distinct magnetic 
ground states, and in particular Neel ordering in spin-1/2 fermions. 
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